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SYMBOLS 


a 


Aa 




C 



F 



m 


n 


^ 1 ^^3 
R 

r 

Ar 

r • X 


T 



column matrix consisting of the coordinates of the position vector 
in the x-*basis 

change in a coordinates between measurements, a 3 x 1 column 
matrix 

range-difference measurements 
orthonormal triplet 

column matrix consisting of a triplet of c vectors 
jk element of 

Grammian of x^, x^, a 2 x 2 matrix 
Grammian of x^^,X 2 , X 3 , a 3 x 3 matrix 
jk element of G~^ 

3x4 matrix consisting of T"^ and row sum of T ^ 

vectors from corresponding reference points to vehicle position 

4 x 1 matrix consisting of product of measurement and the 
corresponding error in the same measurement 

2 x 1 column matrix consisting of the first two elements of 
matrix a 

reference points 

vehicle position 

position vector of vehicle 

change in position vector between different time of measurement 

projection vector; a column matrix consisting of r . x^, r • x^, 

and r • X 3 for simultaneous measurement or of r^^^ • xi, 
f2 j f 3 j 

r^ ^ • x^, r ‘ ^3 nonsimultaneous measurements 

linear transformation between the x and c bases, a 3 x 3 
matrix 

kl element of 


V 



V 

column matrix consisting o£ displacement in the w coordinates 
to account for vehicle motion between measurements 

w 

column matrix consisting of coordinates of the position 

vector 

Aw 

change in w coordinates between different time of measurements 

X 

column matrix consisting of a triplet of x vectors 




position vectors of corresponding reference points 



a 

2x1 column matrix consisting of the first two elements 
projection vector r • x 

of 

the. 

B 

3x1 matrix consisting of coordinates of x^ expressed 
X basis 

in 

the 


coordinates of expressed in the x basis 



0 

angle between line of sight with respect to the horizon; 
numerically positive if measured in the counterclockwise 
direction 


Ar 

longitude of vehicle 



Xr 

latitude of vehicle 



5k 

length of vector 



P 

length of position vector r 




Operators 



c( ) 

cosine of ( ) 



cov( ) 

covariance of ( ) 



d( ) 

derivative of ( ) 



s( ) 

sine of ( ) 



C 

transpose of a matrix 



C ) • ( ) 

inner product of two vectors 



( 

inverse of a matrix 




VI 



Superscripts and Subscripts 


sum of ( ) for all values of i and j, unless specified 
i = 1 , 2 , 3 , and j = 1 , 2 , 3 

index to the time of measurement 

k, Z element of ( ), or k component of ( ) 

See explanation for and 

latitude of vector c^ 



A UNIFIED FORMALISM FOR POSITION DETERMINATION 


FROM VARIOUS NAVIGATION MEASUREMENTS 
Homer Q . Lee 
Ames Research Center 


SUMMARY 


A new computational technique has been developed for coordinate determina- 
tion of moving vehicles. This technique is applicable to a wide class of 
measurements, such as range, range rate, angle, distance sum, distance dif- 
ference, or combinations of these. The generality of the method and the sim- 
plicity of the required computation make the technique ideally suitable for 
navigation of aircraft or spacecraft. 

Coordinate determination by the new technique proceeds in the following 
sequence. First, a vector quantity called the projection is computed by 
simple arithmetic operations on the given measurements. Second, a matrix is 
constructed that depends only on some chosen Cartesian coordinate system and 
on the location of the reference points to which measurements are being made. 
Then the coordinates of the vehicle in the chosen Cartesian system are 
obtained by multiplication of this matrix by the projection vector. The com- 
puted Cartesian coordinates of the vehicle can readily be converted into lon- 
gitude and latitude. A further advantage of this technique is that errors in 
measurements are related directly to errors in computed coordinates by the 
same matrix. 

Different types of measurements affect only the method of constructing 
the projection vector. In this report, construction of the projection vector 
for range measurements, angle measurements, distance-sum measurements, and 
distance-difference measurements is explicitly carried out. It was found that 
the number of measurements, and not the type of measurements , determines the 
uniqueness of the solution. For instance, four measurements, be they range, 
range rate, etc., determine the position uniquely, whereas three measurements 
always yield two solutions to the coordinate-determination problem. In all 
cases, multiple solutions, if present, are exhibited in the process of con- 
structing the projection. 


INTRODUCTION 


The fundamental problem in navigation is determining the position of a 
vehicle by means of measurements relative to known points of reference. Many 
types of measurements can be used. Measurements of range, range rate, angle, 
and differences in distance are commonly employed. Position is determined by 
processing these data in an appropriate manner. 
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The analytical methods used to process the data are as varied as the 
measurements themselves. In one method (refs. 1, 2) position is established 
from the solution of simultaneous quadratic algebraic equations. Another 
method (ref. 3) has the desirable quality of expressing the coordinates of 
position explicitly in terms of measured and known quantities, but requires 
many trigonometric operations that consume a considerable amount of computa- 
tion time. Furthermore, this method is restricted because it assumes that the 
reference points lie in the equatorial plane. Still another method (ref. 4), 
to be used with measurements of distance differences, requires an iterative 
process or spherical trigonometry. This method too is restricted by the 
assumption that the reference points and the vehicle whose position is to be 
determined lie at the same altitude. 

The purpose of this report is to develop another method of processing the 
measurements. This method relies on linear algebra to yield a single formal- 
ism by which any of the usual measurements can be processed but is not limited 
to special navigation situations. Thus, it can be used to determine the posi- 
tion of space vehicles as well as ships and airplanes. It relates the desired 
position parameters directly to the measured quantities through a transforma- 
tion whose terms are already known. 

The idea of the method is very simple and is given in the first section 
below. Any navigation system relies on the use of known points of reference 
for making measurements and on the existence of a standard coordinate system 
for expressing the position of the vehicle. The reference points form a basis 
for a coordinate system and the position of the vehicle is given in this basis 
by the various possible measurements. The reference points are said to be 
known in the sense that their positions are established in terms of the stan- 
dard coordinate basis. Hence the same transformation that relates the basis 
formed by the reference points to the standard also relates the measurements 
to the desired standard expression for vehicle position. The difference 
between the various navigation schemes is shown to reduce to different modes 
of expressing the vehicle’s position in the reference basis. 

The basic scheme works if a set of measurements is made essentially at the 
same time and if the set is sufficient to determine position uniquely. The 
modification to the scheme needed to account for motion of the vehicle between 
measurements is shown to be minor. A more significant modification is then 
developed for the common case where the measurements do not uniquely specify 
the vehicle’s position, but determine it only to within a choice of sign. 

Once this modification is given, the usual navigation systems can be treated. 

Since the formalism developed in this report separates the measurements 
and the transformation between the reference and the standard bases, measure- 
ment errors are similarly separate. The final section of the report develops 
this condition to show how the covariance of position determination and of 
measurement errors are related by a precomputable transformation. 
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DESCRIPTION OF THE BASIC PROCEDURE 


The basic formalism for position determination can most easily be 
developed by reference to the situation illustrated in sketch (a) . Three 

reference points. Pi, P2, and P3 , and 
their positions are located relative to 
an origin 0, say the center of the 
earth, by the vectors xi , X2, and X3. 

The vehicle's position R is located by 
the vector r from 0. The vectors 
I2, and Z3 locate R relative to the 
three reference points. Not illustrated 
are three unit vectors ci , C2> and C3 
forming an orthonormal triplet at 0. 
These vectors form the basis of the stan- 
dard coordinate system; for example, C3 
pointing north from the center of the 
earth, ci lying in the equatorial plane, 
and together with C3, determining the 
Greenwich meridian. 

Since the x's and c's are sets 
of independent vectors, the following 
relations hold: 



t- 

r = a X = w c (1) 

Here the bar over x and c indicate that these quantities are triplets of 
vectors. The a and w are triplets of scalars. The superscript t denotes 
transpose. The a's will be seen to be related to the measurements, and the 
w's are the desired coordinates of the point R. Written out more fully, 
equation (1) appears as 


r = aixi + a2X2 + ^3X3 = wiCi + W2C2 + W3C3 

The two bases, the vector sets x and c, are related by a linear transforma- 
tion 

X = Tc (2) 

The elements in the kth row of T, k = 1,2,3, are the coordinates of x]^ 
expressed in the c-basis. If the reference points are known points on the 
earth, then the elements of T are fixed numbers. If they are moving points, 
such as satellites, then the are known functions of the ephemerides. 

The relation between a and w is 

a^T = ; w = T^a (3) 

Writing equation (1) in the form r = x^a, one obtains 

- _ -t 

r-x=x«xa=Ga ( 4 ) 
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where r • x is the (column) vector (r • Xx , r • X 2 , r . X3) , G is the 


-r 

so-called Grammian x • x , 

and 

can be 

written in 

matrix 

for 


"xi 

• 

X2 • Xi 

X3 • 

x" 

-t 






X • X = 

Xl 

. X2 

X2 • X2 

X3 • 

X2 



• X3 

X 2 • X 3 

X3 • 

X3 


Reference to sketch (a) gives rise to the following equations: 

(r - x^) (r - Xk) = = p2 + 5^2 . 2t ■ ; k = 1, 2, 3 

where = r . r; = xj^ • xj^; and ' ^k* Solving for r • Xj^ 

gives the equation 


r • 




( 5 ) 


In order to illustrate the formalism simply, consider an example in which the 
values of the A^ are given by measurements of range, and in which the alti- 
tude p of the vehicle R is known. Since the are the lengths of the 

vectors to the reference points, they are known. Hence, the right-hand side 
of equation (5) can be evaluated. 


Equations (3) and (4) yield the equations 


w = T^a = T^G ^ (r • x) (6) 

Now the Grammian can be expressed in the following way 

— t — t 

G = XX = Tcc T = TT 
or 

= (TT^)'^ (7) 

_-t 

because cc is equal to the identity matrix. Using equations (6) and (7), 
then, one obtains the equation 

w = f (TT^)'^r . x) = T'^(r • x) (8) 

Equation (8) is the desired result. It shows that the desired coordinates of 
the vehicle w are related to the measurements contained in what will be 
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referred to as the projection vector r • x by a linear transformation whose 
elements are already known. The coordinates of the vehicle can be expressed 
in many ways using the wj^. If its position is to be expressed in terms of 
latitude and longitude for example, and the standard coordinate is 

the Greenwich meridian, then 


X-r = tan 


.1 ^ 
w. 


= tan 


-1 




W<- 


2 2 
1 ^^2 


In any case, the w^ are given in equation (8) by using only arithmetic opera 
tions. Furthermore, the procedure summarized by equation (8), although 
developed with the simplest navigation situation in mind, will be shown to 
hold whenever four independent measurements are given at a time. By the modi- 
fications to the procedure discussed next, any of the usual navigation mea- 
surements can be related to desired coordinate expressions in a unified way. 


MODIFICATIONS TO THE BASIC PROCEDURE 


Nonsimultaneous Measurements 


In deriving the changes to be made to equation (8) when the measurements 
are not given simultaneously, one should keep in mind an example similar to 
that already used. Thus, three fixed reference points are assumed to be given 
and range measurements taken with respect to them. If an airplane is the 
vehicle from which the measurements are made, its motion between range measure 
ments is assumed to be measured. The range measurements are made in the order 
I, 2, then 3. Position will be determined as of the time of the last measure- 
ment. Quantities will be indexed as to the time of measurement by a super- 
scripted number in parentheses. Then the position of the vehicle at the 
various times can be written: 


= r 

- = 

t- 

a X - 


t- 

= w c - 

aw'‘>'8 


= r 

- = 

t- 

a X - 


t- 

= w c - 


> 

(3) 

r^ ^ _ r 

= 

t- 
a X 


t- 
= w c 

> 


As in the previous 

section. 

t- 

a X = 

w^c and 

X = Tc. 

Hence, equations 


(9) 


can be written in the form 
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.(1) _ t 


= a X - Aw 


(Ot- 


( 2 ) t- ( 2)t- 

r =ax-Aw c 


t- -t 
= a X = X a 


( 10 ) 


The projection vector for this case is similar to that given before: 
(r • x)^ = (r^^^ . xi, r*^^^ • X2, r • X3) 


The terms on the right are given by 


.(k) 


Xk = y (P 


2 + £ 2 
k 


- 


wherein the 5^ known, and the and are measured. Then equa- 

tion (10) can be written in the form 


Here 


t 


V 


vi 


(vi, V2, V3) and v^ 

Aw^^^^c . xi 

^Awi^^^Ci + Aw2^^^C2 
Aw^. j = 1, 2, 


-“t 

r * X = XX a - V 
can be written as 

+ Aw 3*-^^C3^ • (tiiCi 

3 


+ ti2C2 + ^1303) 


(11) 


Hereafter, whenever a subscript is repeated it means summation and unless 
specified j =1, 2, 3, that is. 


tj. AW. 


3 

Z 

j=i 


'ij 


Aw. 

J 


The t , , k = 1, 2, 3 are the elements in the first row of the known trans- 

( 1 1 

formation matrix T. The three values refer to the measured motion 

of the airplane between times (1) and [3). The component V2 is similarly 
expressed as 


V 


2 




j = 1, 2, 3 
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The third component V3 is zero. Using equations (6) and (7) one can write 
equation (11) as 

r • X = Tcc^T^(T^)'V - v = Tw - v 
Solving for w, one gets 


w = T ^(r • X + v) (12) 

Equation (12) differs from equation (8) only by the additional term which 
modifies the projection vector by adding the motion of the aircraft between 
the times when range is measured. 


Indirect Expression 

Equation (12) shows that the procedure for expressing the position of the 
vehicle in a desired coordinate system, when the measurements are not made 
simultaneously, is nearly the same as that when all the measurements are made 
at once. The desired coordinates are expressed directly in terms of the mea- 
sured quantities. This direct expression was shown to arise when range mea- 
surements are made to three independent reference points and the altitude of 
the vehicle is known (or measured) . It will be shown later that this form 
holds whenever four independent measurements are made. If only three are 
available, the position of the vehicle can still be determined, but not 
uniquely, and it cannot be expressed quite so directly as before. 

The modification to the basic procedure necessary when only three measure 
ments are available will be developed using the example of two given reference 
points. Range to these points and altitude will be the measurements made. 

Sketch (b) illustrates the situation. 
The reference points and ?2 are 

located by the known vectors x^ and X2 
with respect to point of origin. The 
vehicle is located at R by the vector 
r. The vectors l\ and I2 locate R 
relative to the reference points. The 
lengths of xi and X2, namely and 
are known. The lengths of Zi, l2> 
and r, namely X^, X2, and p, are mea- 
sured simultaneously. 

The two vectors x^ and X2 by 
themselves do not form a basis for spec 
ifying the position of the vehicle. 

They can be used to generate the third 
vector X3 and thus produce a triplet 
X3 is constructed from the cross 



Sketch (b) 

of independent vectors x. The vector 
product of xi and X2' 
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(Xi) X (X2D 


(13) 


X 3 = 


Cxi) X (X2)l 


Since the denominator in equation (13) is the magnitude o£ the numerator, the 
vector X 3 has unit length. 


Equations (1) through (4) hold: 


t- t- 
r = a X = w c 


X = Tc 


(1) 


( 2 ) 


a^T = ; w = T^a 


-_t 

r • X = XX a = Ga 


The two quantities a and n are defined for convenience by 

= (r • xi, r • X 2 ) 

= (ai, a2) 

Then equation (4) can be written 


r • X = 




where F is given by 


F = 


= F 


(3) 

(4) 

(13a) 

(13b) 


(14) 


xi • xi xi • X2 

X2 • Xi X2 • X2 

The following equation comes directly from equation (14) : 

n = F 

Neither r • X 3 nor a 3 is known directly. Either can be calculated by 
considering r . r: 


(15) 


(16) 


r • r = = a 


'xx‘a = [„b a3][^ 


(17) 
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Hence, follow the equations 


a3^ = - n^Fn 

= - a^F 

or 

a-3 = ±[P^ - (18) 

Since the correct sign must be chosen for in equation (18) , the solution 

is not unique. The choice of sign means that R could be at either of two 
places located symmetrically with respect to the plane formed by the vectors 
X 2 and X 2 * 

Equation (18) completes the information required for position determina- 
tion. The steps are, in summary: 


n = F"^a 

(16) 

a3 = ±[p^ - a^F 

(18) 

w = T^a 

(3) 


Whereas equation (8) expresses w directly in terms of the measurements in 
r • X, equations (16), (18), and (3) form an indirect expression, requiring 
calculations of intermediate quantities. 

If T • xi and r • X 2 are measured at different times, equation (16) 
changes to the form 

n = F + V 

where 

Vi = (j = 1, 2, 3) 

and 

V2 = 0 


CONSTRUCTION OF PROJECTION VECTORS FOR VARIOUS NAVIGATION SYSTEMS 


An important advantage of this formalism is that different types of 
measurements affect only the method of constructing the projection vector, a 
component of which is expressed in equation (5) . Once the equations that 
relate the projection vector and the measurements are known, the basic proce- 
dure or its modification can be applied to solve all position determination 
problems. In the previous section, p and X were the measured quantities in 
the projection. This section is concerned with constructing the projection 
when other types of measurements are used. 
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The following list classifies navigation systems according to the type of 
measurements provided to construct the projection: 

1. Altitude, and any combination of two or three range or angle 
measurements . 

2. Four independent range measurements. 

3. One altitude and three distance difference measurements. 

4. One altitude and two distance difference measurements. 

5. Three independent range measurements, where the plane containing the 
reference points either (a] does not contain the origin, or (b) does contain 
the origin. 

Calculating the projection vector from the data of each of these systems 
will be shown in the order listed. 


Projection Using Angle Measurements 


Equation (5) expresses the projection vector in terms of the known 
altitude p and ^ and the range, A. It will now be shown how an angle measure- 
ment can be substituted for range. Sketch (cl) shows that the equation 

r - Z = X holds. Taking the dot product of x 
with itself and rearranging terms give the 
equations 



0 = - 2 pXc( 90 ° + 0) + p2 - ^2 

= x2 + 2pXs0 + p2 - 52 


(19a) 


where 0, the angle between the horizon and the 
line of sight, is considered positive if the 
reference point, P, lies above the horizon. Solv- 
ing equation (19a) for X gives 


X = -ps0 ±[(PS0)2 + (§2 _ p2jj 


1/2 


(19b) 


Note that A, being a length, must always be posi- 
tive. To choose the correct sign in equation (19b) 
one must consider the two cases 5 > P and 5 < p. Considering first the con- 
dition 5 > p, it will be shown that 


-ps0 + [(ps0)2 + (g^ 




(19c) 
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both for 6^0 and for 0^0. It can be seen from sketch (c2) that the 
following relations hold: 



In both cases where 0 ^ 0 or 0 ^ 0, equation (20) is the required solution 
for X for the condition ^ > p. If, on the other hand, C < P, then there 
are two possible solutions for X. Sketch (c3) and equation (20) show that 
the following relations hold: 


QR = |ps0| = -ps0 since 0 < 0 

PR = ^ - PQ = -pse - /(ps0)2 + ($2 - p2) 

= ^ + pT^ = -PS0 + /(ps0)2 + (52 _ p2) 

Therefore, there are two possible solutions: 

X = _PS0 + [(PS0)2 + (52 _ p2)]l/2 

or 

X = -ps0 - [(ps0)2 + (^2 _ p2)]l/2 

To know which solution to choose under the 
condition ^ < p requires additional 
information. 


R 



Sketch (c3) 


To calculate the projection vector in terms of angle measurements one can 
”dot" through the equation r - t = x with r to get r • x = p^ - pXs0. Then 
using equation (19b) gives 


r • X = p^c0 ±ps0(C^ - p^c^0)^^^ (21) 

After the proper sign in the equation is selected, equation (21) is a variation 
in form of equation (5). Specifically, the equation that relates the kth 
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component of the projection vector and its corresponding angle measurement is 
given by 


(r . 


'p2c0j^ + 

,p2c0j^ - pS0j^[5^ - p2c0j^]^/2 


(5a) 


Projection Using Four Independent Range Measurements 

When measurements of range are given with respect to four independent 
reference points, the following procedure can be used to calculate the projec- 
tion vector. Let , X 2 , X 3 , and X 4 be the position vector of the four 
reference points relative to a common origin. One of them, say X 4 , can be 
written as a linear combination of the other three; 


X4 = e^x^ + ^3^3 ~ 

Since all the x’s are known, the g*s are also known. The problem is to 
eliminate the unknown altitude, p, from the expression for the projection. 
Now equation (1) relates the vehicle’s position to three reference points by 
the expression r = a"^x. The projection onto X 4 then is given by 


r . x^ = a"^xx^B = (x • x )^3 

where equation (4) has been used. This equation can be expanded by means of 
equation (5) to give 


— ^ 

t 


- A, 2 , 


3i 

p2 + ^^2 - X2^ 


32 

r.2 . r 2 ), 2 

_P * ^3 " ^3 


^3 


The solution of this equation for is 


C23a) 






- 




(23b) 


The denominator of equation (23b) cannot be zero so long as the four reference 
points do not lie on a plane. After p^ is determined from equation (23b) 
all the quantities in the right side of equation (5) are known and the projec- 
tion vector is completely specified. 
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Projection Using Altitude and Three Distance Difference Measurements 


Distance-difference measurements arise in using hyperbolic intersection 
methods of navigation such as Loran C. Sketch (d) shows the point R to be 

determined and four reference stations through 
P 4 . These known points are at positions 
through X 4 from an origin. A constant difference 
in distance of R from two reference stations 
gives the locus of a hyperbola. This difference 
is determined by measuring the time difference 
between pulses from each station. The length ^ 
of the vector to the first reference of any pair of 
stations is not known. Only the difference in 
their distances is known. With Pi as the primary 
references, sketch (d) illustrates the following 
relations between the vehicle at R and four 
known stations: 



= X ' 


X^ - X^ + b^ 


A3 - Ai + 

A^^ = Ai + b3 


(24) 


As in equation (22) of the previous section, x^ can be written in terms of 
the other vectors 


= 3^X3 + 63X3 + 33X3 

where the 6 ’s are known. Equation (23a) results. Substituting for the Xj^^s 
by equation (24) and expanding gives the following quadratic equation for X: 

0 2X(32b^ + - b^) +^3j_(p^ + ?^^) - (p^ + 5^^) 

+ b^^ + b^^ + b 3 ^ = 0 (25) 

After X is computed from equation (25) , then equation (24) can be used to 
compute Xi, X 2 , and X 3 . Therefore, all quantities in the right side of 
equation (5) are known and so is the projection vector. 


Projection Using Altitude and Two Distance Differences 

If one altitude and two distance differences are known, the equation for 
X becomes a quartic. Two solutions are possible; the proper one must be 
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selected. The equation to be solved is 



p2 + 5^2 _ ^^2 


p2 + 5^2 _ ^^2 

>— 1 

II 

CM 

Q. 

p2 + ^2^ - + bj)2 

G'^ 

p2 + - (Aj + b^)2 


_P^ + €3^ - (Aj + b2)^. 


^P^ + ^ 3 ^ - (Ai + ^2)^. 


( 26 ) 


Equation (26) comes from equations (1) and (4) and taking a dot product 


r = a^x . x\ = a^Ga = (r • x)^(G ^^GG ^ (r • x) 


(r • x)^G ^(r . X) 


With the measured quantities p, b^, and hz, and the notation 
the ijth element of G and defining the matrix P as 


g . \ for 
ij 


P = 


p 2 + 532 . b 2 ^ 


equation (26) can be written as the quartic 


^‘'2 g-j - 4x3[b3,b2] 


IJ 


+ 4X2/[bj,b2] 


-1 

1 — 1 

1 

-r 

§21 

^22 

®23 

-1 

-1 

-1 

_^31 

®32 

^3 3_ 


-1 

-1 

®22 

^23 

-1 

-1 

^32 

^33 


’ 2^2i ’ ^^3ij 


> P 


4A[b^,b2l 


““ ”1 

-1 

-1 

§21 

®22 

®2 3 

_1 

-1 

-1 

^31 

®32 

^33 


P + P^G ^P - 4p‘ 


0 (.27) 


Computation of the coefficients for equation (27), although tedious by- 
hand, requires a simple digital computer program involving only matrix 
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multiplication and basic arithmetic operations. After the coefficients are 
computed, methods, such as Ferraries or Descarte*s are available for extract- 
ing the required root for equation (27). Again after X is computed from 
equation (27), the first three expressions in equation (24) can be applied to 
compute Xx, X 2 > and X 3 . This computation completely specifies all quantities 
in the right side of equation (5) . 


Projection Using Three Independent Range Measurements 

Two cases will be considered here. In the first, the reference points 
form a plane that does not contain the origin. The three range measurements 
Xi, X 25 and X 3 are given, and p is to be determined. Then equation (26) 
can be expanded to give the following biquadratic in p : 


p'x:* 

ij 


ij 


+ 2 


{Ehi 



+ 2 





(28) 


In the second case, three range measurements are made again and p is to 
be determined, but this time the plane formed by the reference points does con- 
tain the origin. This development is similar to that in the section labeled 
"Indirect Expression" wherein the reference vectors were not linearly indepen- 
dent. In the present case, only two of the reference vectors are independent. 
Expressing the third in terms of the first two gives the expression 

X3 = (3x^32) (X3^,x^)^ 

Following the arguments used in deriving equation (23a) gives the following 
expression for the projection of r along X 3 : 

[(p^ + ?x^ “ > (P^ + " ^2^)] [ 3 x > 32 ^^ = P^ + $3^ ' X32 


Solving this equation for p' 


’2 
gives 


= [SeiCq" - - (53" - ^3") /(i -S^i) 


i = 1, 2 


(29) 
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Since the reference points are not colinear, the denominator of equation (29) 
is not zero. Again, equation (29) can be used to compute p^. Once is 

knovm, the right side of equation (5) and thus the projection vector is 
completely specified. 


Table 1 summarizes all typical coordinate-determination problems and 
their corresponding solution. 

TABLE 1.- TYPICAL COORDINATE-DETERMINATION PROBLEMS 


' Given set of measurements 

Number of 
solutions 

Brief description of 
solution 

One altitude measurement and 

(a) 3 independent ranges 

(b) 2 independent ranges and 1 angle 

(c) 1 independent range and 2 angles 

(d) 3 independent angle measurements 

1 

Direct expression 
Use equations (5) and 
(21) to construct the 
projection vector 

Four independent range measurements or 

(a) 3 independent ranges and 1 angle 

(b) 2 independent ranges and 2 angles 

(c) 1 independent range and 3 angles 

1 

Direct expression 
Use equation (23b) to 
solve for and use 

equations (21) and (5) 
to solve variations of 
this problem 

One altitude and 

(a) 3 distance differences 

(b) 3 distance sums 

(c) 2 distance sums, 1 distance 

difference 

(d) 1 distance sum, 2 distance 

differences 

Also other combinations of distance 
difference, distance sum, and range 
rate 

1 

Direct expression 
Use equation (25) to 
solve for 

Three independent ranges 

(a) 2 ranges, 1 angle 

(b) 1 range, 2 angles 

2 

Direct expression 
Use equation (28) to 
solve for p^, and 
equation (21) for 
variation of this 
problem 

One altitude and 

(a) 2 distance differences 

(b) 1 distance difference, 1 distance 

sum 

(c) 2 distance sums 

Also other combinations of distance sum, 
distance difference, and range rate 

2 

Direct expression 
Use equation (27) to 
solve for X^ 

One altitude and 

(a) 2 ranges 

(b) 1 range, 1 angle 

(c) 2 angles 

2 

Indirect expression 
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RELATION BETWEEN COVARIANCES OF NAVIGATION AND MEASUREMENT ERRORS 


The procedure developed in the preceding sections converts measured 
quantities into the quantities desired for navigation by means of a transfor- 
mation whose terms are already known. Since the procedure separates the mea- 
surements and the transformation, it also separates measurement errors from 
the particular geometry summarized by the transformation. The separation will 
be shown in this section wherein expressions for the navigation errors and 
their covariances will be related to measurement errors and their covariances. 


Error Covariance Matrix 

The objective here is to show that error of the computed coordinates is a 
linear transformation of the errors in the measured quantities. Furthermore, 
the error covariance matrix of the computed coordinates is also some linear 
combination of the error covariance of the measured quantities. The relation 
between errors in computed coordinates and errors in measured quantities for 
the one altitude and three range measurements case will be derived. Ass\jming 
that the position of the reference points is completely known, we may take 
the partial derivatives of the computed coordinates in equation (8) with 
respect to the measured quantities as follows: 






dA^ 

dw = 


p dp - T ^ 

A2 

dA2 




_^3 

dAg 


Next, we define the following two matrices 


(30a) 


J = 





p-l 


(30b) 


m (p dp ,— AjdAj^, A^dA^,— A^dA^) (30c) 

Equation (30a) may be rewritten in terms of the two matrices defined in equa- 
tions (30b) and (30c) as 

dw = Jm (30d) 
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The error of computed coordinates, therefore, is a linear transformation of m, 
the product of the measured quantities and the error in measurement of the 
corresponding quantity. The linear transformation matrix J depends only on 
the location of the reference points . By virtue of equation (30d) , the 
covariance matrix of the computed coordinates is 

cov(dw) = J cov(m3J^ (30e) 

In many coordinate determination problems, the desired coordinates are lati- 
tude and longitude. Then the rectangular coordinates must be converted into 
polar coordinates. For the direct expression, if the c-vectors are chosen 
such that c^, C 3 are vectors from the center of the earth to Greenwich and 
to the North Pole, respectively, then the longitude and latitude of the air- 
craft position are, respectively. 


-1 


A ^ -1 '*'3 

A = tan — — — — — 


(31) 


The error in computed latitude and longitude is 


dA 


dA 


cos A 


A 


^ 


-w- sin A 
3 r 


-sin A 


/ 




cos X 
3 r 


A 


2 -H 


dw 


(32) 


Again, we define the matrix L 


L = 


cos X 


A 


2 + W 2 


-W3 sin X 


-sin X 


AAA 


-Wo COS X 
6 Y 


/w^^ + W^^ 


(33) 


Next, we write equation (32) in terms of equations (30b) and (33) 


dA 


dA 

u r-J 


= LJm 


(34) 


18 



I 


The error covariance matrix for the computed latitude and longitude are 
related to the error covariance matrix of the measured quantities in the fol- 
lowing expression 


= LJ cov(m)J^L^ (35) 


It is to be emphasized that matrix J in equations (30d) , (30e) , (34) , and 
(35) depends only on the location of the reference points and that the matrix 
L in equations (34) and (35) , which becomes singular at both the north and 
south poles, depends on aircraft position. 



The relation between the errors of computed coordinates and the errors of 
measured quantities will be derived for the indirect expression. First, write 
equation (3) in component form in light of equations , (16) , (13a) , and (5) and 
as follows: 






"r • 






^2 

= T 

0 


r • X 2 

-^3- 


. 0 0 1^ 


- - 



F 


o 

p2 


' ^1.^' 

ft 



0 

p2 


“ ^2^ 


0 

0 

i_ 


^3 



(36) 


Second, take the partial derivatives of the computed Cartesian coordinate in 
equation (36) with respect to the measured quantities 



j 



P dp 

dw = T^ 




-Xg dXj 


1 - a^ - a^ -aj 

-^2 


— X 2 dX 2 


ag ag 

^3 



(37) 


Next we define the matrices 
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-1 

ij 


K = 


j 2J 


1 - a2 - 


F"^ 



(38a) 


m' = (p dp, -Xj dXj, -X 2 dX 2 ) (38b) 

Then upon substituting equations (38a and b) into equation (37) , we have 

dw = T^Km' (39) 

As exhibited by equation (39) , the error of the computed coordinates is 
again a linear transformation of errors in measured quantities. However, in 
this indirect expression case, the linear transformation matrix depends on 
both the location of the reference points and on the position of the aircraft. 
The covariance matrix of the computed coordinates is 

cov(dw) = T^K cov(m')K^T 

To express the computed coordinates for the indirect case in latitude and 
longitude, choose c^ as the vector pointing from the center of the earth to 
the intersection of the orbital plane with the equatorial plane so that the 
longitudinal distance from Greenwich and the projection of c^ on the surface 
of the Earth is no greater than 180°. Choose C 3 as in the direct expression 
case. Then the latitude and longitude of the aircraft are, respectively. 


-1 "2 


A = tan 
r 


-1 


w. 


A 


^ + W2^' 


The errors in computed longitude and latitude are related to the errors in 
measurements as 

rx. 


L-^r- 


= LT‘'Km* 


The covariance matrix of the computed coordinates is given by 

'dX 


cov 




LT^K cov(m')K^TL^ 
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CONCLUSION 




The coordinates of the position of a vehicle, say an aircraft, can be 
determined by a matrix multiplication of the projection vectors constructed 
from measurements. Measurements may be range, range rate, angle, distance 
difference, etc.^ and methods for construction of the projection vector from 
typical combinations of measurements were presented. The matrix displaying 
the coordinates by a linear transformation of the projection vector depends 
on the reference points from which the measurements were made. Since the 
position of the reference points is already known, this matrix is 
precomputable. 

The basic procedure developed for coordinate determination assumes that 
measurements were made essentially at the same time and that the set of mea- 
surements is sufficient to determine the vehicle position uniquely. The 
scheme also works for nonsimultaneous measurements, provided a component is 
added in the projection vector to account for motions of the vehicle between 
measurements. If the set of measurements is insufficient to determine the 
position uniquely, then the procedure can determine the position to within a 
choice of sign. 

The advantages of this technique are: (1) it is applicable to any com- 

bination of different types of measurements; (2) the measurements and the 
transformation are separated and related by a linear transformation; (3) the 
error in computed coordinates and that in measured quantities are separated 
and related by another linear transformation; (4) the covariance matrix of 
the computed coordinates and that of the measurements are also separated and 
related by some variation of the linear transformation in (3); and (5) the 
transformation in (2) is already known and is therefore precomputable. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 9, 1971 
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